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Abstract 



The Riemann-Roch theorem on a graph G is related to Alexander duality in com- 
binatorial commutive algebra. We study the lattice ideal given by chip firing on 
G and the initial ideal whose standard monomials are the G-parking functions. 
When G is a saturated graph, these ideals are generic and the Scarf complex is 
a minimal free resolution. Otherwise, syzygies are obtained by degeneration. We 
also develop a self-contained Riemann-Roch theory for artinian monomial ideals. 

1 Introduction 

We examine the Riemann-Roch theorem on a finite graph G, due to Baker and Norine 
[3], through the lens of combinatorial commutative algebra. Throughout this paper, 
G is undirected and connected, has n nodes, and multiple edges are allowed, but we 
do not allow loops. Its Laplacian is a symmetric n x n-matrix with non-positive 
integer entries off the diagonal and kernel spanned by e = (1, 1, ... , f). Divisors on G 
are identified with Laurent monomials x u = The chip firing moves are 

binomials x u — x v where u, v > and u — v is in the lattice spanned by the columns of 
Aq- The lattice ideal Iq spanned by such binomials is here called the toppling ideal of 
the graph G. It was introduced by Perkinson, Perlman and Wilmes [TH [15], following 
an earlier study of the inhomogeneous version of Iq by Cori, Rossin and Salvy [6]. 

For any fixed node, the toppling ideal Iq has a distinguished initial monomial ideal 
Mq. This monomial ideal was studied by Postnikov and Shapiro [12] . and it is character- 
ized by the property that the standard monomials of Mq are the G-parking functions. We 
construct free resolutions for both Iq and Mq, and we study their role for Riemann-Roch 
theory on G. For an illustration, consider the complete graph on four nodes, G = K4. 
The chip firing moves on K± are the integer linear combinations of the columns of 



The toppling ideal is the lattice ideal in K[x] = K[xx, 22, ^3, £4] that represents image z (Ac) 
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This ideal is generic in the sense of Peeva and Sturmfels [TU], as each of the seven binomi- 
als contains all four variables. The minimal free resolution is given by the Scarf complex 

< — K[x] < — K[x] 7 < — K[x] 12 < — K[x] 6 < — 0. (3) 

The seven binomials in ^ form a Grobner basis of Iq, with the underlined monomials 
generating the initial ideal Mq- That monomial ideal has the irreducible decomposition 

M G = (x t ,xl, xl) n (xt,xl, xl) n (xl, x 2 , xl) n (xj, x\, x 3 ) n (xl, x 2 , x\) n (xl, x\, x 3 ). 

The ideal Mq is the tree ideal of [U §4.3.4]. Its standard monomials are in bijection with 
the 16 spanning trees. Its Alexander dual is generated by the six socle elements 

2 2 2 2 2 2 ( a\ 

*^2*^3; *^2*^3? *^1*^35 *^1*^2' ^l^S? X-^X 2 . v / 

These correspond to the maximal parking functions studied in combinatorics; see |5],I12|. 
We claim that the duality seen in Figures 4.2 and 4.3 of [9] is the same as that expressed 
in the Riemann-Roch Theorem for G. This will be made precise in Sections 3 and 4. 

The present article is organized as follows: Section 2 is concerned with the case when 
G is a saturated graph, meaning that any two nodes i and j are connected by at least one 
edge. We show that here Iq is a generic lattice ideal, and we determine its minimal free 
resolution and its Hilbert series in the finest grading. The Scarf complex of the initial 
monomial ideal Mq is supported on the barycentric subdivision of the in — 2)-simplex 
[j~2| §6], and this lifts to the Scarf complex of the lattice ideal Iq by pi), Corollary 5.5]. 

In Section 3 we revisit the Riemann-Roch formula 

rank(.D) — rank(fT— D) = degree ( D )— genus + 1. (5) 

We prove this formula in an entirely new setting: the role of the curve is played by a 
monomial ideal, and that of the divisors D and K is played by monomials x b and x K . 
The identity (|5]) is shown for monomial ideals that are artinian, level, and reflection- 
invariant. This includes the parking function ideals Mq derived from saturated graphs G. 

In Section 4 we extend our results to the case of graphs G that are not saturated, 
and we rederive Riemann-Roch for graphs as a corollary. Here Mq is still an initial 
ideal of Iq, but the choice of term order is more delicate [HI §5]. One choice is the cost 
function used by Baker and Shokrieh for the integer program in [U Theorem 4.1]. The 
Scarf complexes in Section 2 support cellular free resolutions of Iq and Mq, but these 
resolutions are usually far from minimal. We conclude with several open questions. 

This paper demonstrates how Riemann-Roch theory embeds into combinatorial com- 
mutative algebra. Our main results are Theorems [2| 13 and 25 These build on earlier 



works, notably [TJ and [12], but they go much further and are new in their current form. 

When this collaboration started in the summer of 2011, both authors were unaware of 
the articles [H1IT5] written on similar topics by David Perkinson and his students at Reed 
College. As our point of departure, we chose to focus on chip firing in the most classical 
case of undirected graphs, but with the tacit understanding that our ideals and modules 
generalize to directed graphs, arithmetic graphs, simplicial complexes, matroids, abelian 
networks, or any of the other extensions seen in the recent chip firing literature (cf. [2]). 
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2 Saturated graphs 



In this section, we assume that the graph G has Uy edges between node i and node 
j, where is a positive integer, for i ^ j. However, we do not allow loops, so that 
!iu = U22 = ■ ■ ■ = u nn = 0. Thus, in the language of [T2"| . G is a saturated graph. 
We shall see that, under this hypothesis, the lattice ideal Iq is generic, and an explicit 
combinatorial description of its minimal free resolution can be given. Throughout this 
paper we work in the polynomial ring K[x] = K[xi, . . . , x n ] over an arbitrary field K. 

We begin by explicitly showing the generators of the lattice ideal Iq in the case n = 4. 
Example 1. If G is a saturated graph on [4] = {1, 2, 3,4} then Iq is generated by 

u12+u13+m14 ™m12™m13™«14 ™«12+m23+w24 ™«12 ™m23 ™«24 ™«13+'"23+m34 ^,"13 ™«23 „u34 

Jj-j^ "^2 3 4 5 2 1 3 4 5 3 1 2 4. 5 

™«14™M24™«34 ™«14+M24+«34 „,U13+«14 „,«23+«24 ™«13+«23 ™«14+M24 

1 2 3 "^4 5 1 2 3 *^4 5 

W12+W14™^23+W34 _ ^^12+^23^^14+^34 ^12+^24^13+^34 _ ^12+^13 U24+W34 

Here the i/y are arbitrary positive integers. These binomials form a Grobner basis. The 
initial ideal Mq is generated by the underlined monomials. The minimal free resolution 
of Iq has the form ([3]). The same holds for Mq, as was shown in [T21 Corollary 6.9]. The 
minimal resolution of Mq is given by the Scarf complex, which is depicted in Figure[T} □ 

We now state our main result in this section. For disjoint subsets I and J of [n] we set 

n 



A split of the set [n] = {1,2, ... ,71} is an unordered pair (/, J) of non-empty disjoint 
subsets / and J whose union equals [n]. The number of splits equals 2 n_1 — 1. With 
each split (/, J) we associate the following binomial which is well-defined up to sign: 



These are precisely the seven binomials in Example [TJ one for each split (/, J). 

Let Cyc n k denote the set of cyclically ordered partitions of the set [n] into k blocks. 
Each element of Cyc n k has the form (ii, I2, . . . , Ik), where I\ U I2 U • • • U Ik = [n] is a 
partition. We regard the (ii, h, ■ ■ ■ , h) as formal symbols, subject to the identifications 

(ii, J2, . . . , Ik-i, Ik) — (h, h, ■ ■ ■ , h, h) = ••• = (h, h, ■ ■ ■ , h-2, h-i)- 

We write K[x] Cyc ™> fe for the free K[x]-module generated by these symbols. The rank of 
this free module equals the number of cyclically ordered partitions, namely 

|Cyc njfc | = (k-l)\-S n>k , (7) 

where S n ^ is the Stirling number of the second kind, i.e., the number of partitions of 
the set [n] into k blocks. Let CyCc denote the following complex of free K[x]-modules: 

^~ KM Cyc '^ <— K[xl Cyc -- 2 < — K[xl Cyc - 3 < — ■■■ < — K[xl Cyc — < — 0, (8) 
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(i^+u^+^.O.O) ("13 + ^4, "23+ "24'°) ( > ^2 +u 23 +u 24 > 0) 

Figure 1: The barycentric subdivision of the (n — 2)-simplex supports minimal free 
resolutions for toppling ideals of saturated graphs with n nodes and their initial ideals. 
Shown here is the Scarf complex d3]) for n = 4, with 7 vertices, 12 edges and 6 triangles. 



where the boundary map from K[x] Cyc ™> r to Rfx] 070 ™' 1 — 1 is given by the formula 

(Jl ' /2 ' J3 '---' Jr) * -x^.(/ 2) / 3 ,...,/ r -i,/iUJ r ). • (9) 

In this formula it is assumed that n G 7 r , so as to ensure that all signs are consistent. 

Theorem 2. Lei G be a saturated graph. The toppling ideal Iq is a generic lattice ideal. 
It is minimally generated by the 2 n ~ 1 — 1 binomials |^) ; these form a reverse lexicographic 
Grobner basis, the complex CyCc coincides with the Scarf complex, and this complex 
minimally resolves K{x\/Iq. 

Proof. We begin by noting that x /_5-J — x J ^ 7 actually lies in the ideal Iq. To see this, 
let ei denote the incidence vector in {0, l} n that represents the subset I of [n\. The z-th 
coordinate of the vector Aq ■ ej is equal to J2ke j u ^ it i €. I, and it is — J2kei Uik ^ 
i G J. Hence Ac ■ ei is represented algebraically by x 7 ^ J — x J ^ / , which is hence in Iq. 

Fix any reverse lexicographic term order on K[x] that has x n as the smallest variable, 
and let in(Jfj) denote the initial monomial ideal of Iq- Since Iq is a lattice ideal, x n 
is a non-zerodivisor and it does not divide any of the generators of in(/(j). We may 
thus regard in(/(j) as an artinian ideal in K[x\ n ] = K[xi, . . . ,x n _i]. The index of the 
Laplacian lattice image z (Ac) in its saturation {u G Z n : u\ + • • • + u n = 0} equals the 
number Tq of spanning trees of G. Hence in(/c) has Tq standard monomials in K[x\ n ]. 

Let Mq denote the ideal generated by the initial monomials of the binomials in (|6]): 

Mq = ( x 7 -+M\-r : J non-empty subset of [n - 1] ). (10) 

By construction, the inclusion Mq C iti^Iq) holds. The monomial ideal Mq was studied 
in [12] and shown to have precisely Tq standard monomials. Indeed, the standard 
monomials of Mq are in bijection with the n-reduced divisors. It is known in the chip 
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firing literature (cf. [31 OS IB]) that their number equals the number Tq of spanning trees. 
Hence Mq and in(Ic) are artinian of the same colength in K[x\ n ], so they must be equal: 

M G = in(/ G ). 

Therefore the binomials (|6| form a Grobner basis, and hence a generating set, of Iq. 

The ideal Iq is a generic lattice ideal, in the sense of [TO], because all n variables 
xi,...,x n occur in the binomial ([6]). Here we are using that G is saturated. By [TUl 
Theorem 4.2], the Scarf complex is the (essentially unique) minimal free resolution of Iq. 

It remains to be seen that the Scarf complex is equal to CyCc- Postnikov and 
Shapiro [121 Corollary 6.9] showed that the Scarf complex of the initial ideal Mq is 
supported on the barycentric subdivision of the (n — 2)-simplex, as shown in Figure [I] 
The Scarf resolution has the format (Jsp , but with K[x] replaced by K[x\ n ]. Here, we 
label the cells in that barycentric subdivision with ordered partitions (Jx, J 2 , . . . , J r ) 
satisfying n G I r . The boundary maps in the Scarf resolution are then given by Q, 
namely, by the sum ranging from s = 1 to s — r — 1, but without the additional term 

-X^ 71 ■ (/ 2 ,/ 3 ,..., JiU/r). 

We pass from the Scarf resolution of Mq to that of Iq by the combinatorial rule in 
[TUt Theorem 5.4]. This adds precisely one term to the boundary of each Scarf simplex 
of Mq. In our case, that additional term is precisely the one above, and we get (|9]). □ 

Example 3. Returning to Example [TJ with the seven binomials in that order, here are 
the matrices over K[xi, x 2 , x 3 , x 4 ] that represent the first and second syzygies in CyC^. 





(1,2,34) 


(2,1,34) 


(1,3,24) 


(3,1,24) 


(2,3,14) 


(3,2,14) 


(1,23,4) 


(12,3,4) 


(1,234) 


/_ x 2^34 


_ x 34->-2 


_ x 3^24 


_ x 24^3 








_ x 23^4 .. 
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(2,134) 


_ x 34^1 


_ x 1^34 








_ x 3^14 


_ x 14^3 








(3,124) 








_ x 24^1 


_ x 1^24 


_ x 14^2 


_ x 2^14 





■ -x 4 ^ 12 


(123,4) 




















x™ • • 


• x 12 ^ 3 


(12,34) 


x 1^2 


x 2^1 

















• -x 3 ^ 4 


(13,24) 








x 1^3 


x 3^1 














(23,14) 


V o 
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x 3^2 


-x 4 ^ 1 • • 


• o / 





(1,2,3,4) 


(1,3,2,4) 


(2,1,3,4) 


(2,3,1,4) 


(3,1.2,4) 


(3,2,1,4) 


(1,2,34) 
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Note that the seven binomial generators of Iq appear as the 2 x 2-minors of the 3x2- 
matrices seen in the six pairs of columns within the 7 x 12-matrix of first syzygies. The 
syzygies of the ideal Mq generated by the underlined monomials in Example [T] are found 
by replacing with all monomials that have the symbol "4" to the left of the arrow. □ 

One immediate application of our minimal free resolution is a formula for the Hilbert 
series of the ring K[x]/Jg in its natural grading by the group Div(G) = Z^/image^Ac). 
As is customary in chip firing theory [2J [3j HJ E] , we consider the decomposition 

Div(G) = Z © Divo(G), 

where Z records the degree of a divisor on G, and Div (G) is the finite subgroup of 
divisors of degree 0. The order of Divo(G) is the number of spanning trees of G. Let 
t and q denote the generators of the group algebra Z[Div(G)] corresponding to this 
decomposition. The Hilbert series of K[x]/Jg equals 1/(1 — t) times the Hilbert series 
of K[x\ n ]/M<3, where M G = in(J G ) is the initial ideal in (10). The latter series equals 



V^ t |u| ? div(u)_ 
u 

This finite sum is over all elements u e N n_1 that represent parking functions on G with 
respect to the last node n, and div(u) denotes the class of the reduced divisor of degree 
given by the vector (u, — ^2ui). See [TH Theorem 6.14] for a nice formula, due to 
Merino [8] , which expresses this sum with q = 1 in terms of the Tutte polynomial of G. 

We fix the natural epimorphism from the semigroup algebra of to that of Div(G): 
i) : Z[x\ n ] -> Z[Div(G)] , x u h-> t^q^ u \ 
With this notation, our minimal free resolution in Theorem [2] implies the following result: 
Corollary 4. The Hilbert series of K\x\/ Iq in the grading by the group Div(G) equals 

1 \- Au\ div(«) 1 ~ ELll- 1 ^ E(J 1 ,J 2 ,...,J fc )eCyc w , fc VKx Jl ^ /2 X J ^ J3 ■ ■ ■ X /fc -^) 

q " {l-t){l-^{x x )){l-^{x 2 )) ■■■ (l-V(Zn-l)) 

Proof. It suffices to note that the Z n_1 -degree of the basis element (Ji, I 2 , . . . , h) of the 
free K[x\ n ]-module in the k-th step of the resolution of Mq is the exponent vector of 
x Jl ^ /2 x /2 ^ /3 • • .x /fe - 1_i>/fe . This monomial does not contain x n because n G Ik- □ 

We close this section with a combinatorial recipe for the socle monomials modulo Mq. 
These are the monomials x u that are not in Mq but x 11 ^^ e Mq for i = 1, . . . ,n — 1 . Each 
permutation of [n — 1] corresponds to a flag T of subsets C T\ C T 2 C T 3 C • • • C T n _ x . 
The flag is complete, meaning that each inclusion is strict and each Tj\Tj_i is a singleton. 
Let Tj denote the set complement of T, with respect to [n]. For instance, T n _i = {n}. 

Corollary 5. T7ie soc/e monomials o/K[x\„]/Mq are precisely the (n — 1)! monomials 

S T = lcmfx^^ 1 , X T2 ^ 2 , . . , X^-^^-'j/fxHs • • • (11) 

where T runs over all complete flags of subsets of [n — 1] . 
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Proof. The Scarf complex of Mq is a minimal free resolution and it is supported on the 
barycentric subdivision of the (n — 2)-simplex, by (T2J Corollary 6.9] and our discussion 
above. Each facet in that barycentric subdivision corresponds to a complete flag T. The 
vertices of that facet are labeled by x Tl_>Tl , x T2->T2 , . . . , x 7 ™- 1-5 * 7 ™- 1 in the Scarf complex, 
and the monomial label of the facet is their least common multiple. Facets of the Scarf 
complex are in bijection with the irreducible components of Mq and also with the socle 
monomials modulo Mq. By [9, Corollary 6.20], each socle monomial is multiplied by the 
product of all variables to give the monomial label of the corresponding facet. □ 

Remark 6. Our results hold verbatim for all generic sublattices of finite index in the root 
lattice A n = {u G Z n : Y^7=i Ui = so we recover the Voronoi theory of [HE]- We posit 
that our commutative algebra derivation of their Voronoi theory is a natural and useful 
one, and that it opens up new and unexpected connections. For instance, Grobner bases 
of lattice ideals are fundamental for integer programming [13]. One original source for 
that application is Herbert Scarf's seminal work on neighborhood systems in economics. 
A key example that motivated Scarf was the Leontief system [T31 §2A]. It turns out 
that the lattices representing Leontief systems are precisely our generic lattices here. 
The Grobner basis property stated in Theorem [2] is in fact equivalent to [131 Theorem 
2.2]. ~ □ 



3 A Riemann-Roch Theorem for Monomial Ideals 

In this section we fix an arbitrary artinian monomial ideal M in a polynomial ring 
K[x] = K[xi, . . . ,x m ). We focus on Alexander duality [9J §5], and we establish the 
Riemann-Roch formula ^ in this new context. Towards the end of this section, and 
in the next section, we will recover the Riemann-Roch formula for graphs from the 
Riemann-Roch formula for monomial ideals. To begin with, we need to gather the 
ingredients, that is, we need to redefine the notions of divisor, genus, rank and degree. 

The role of divisors on the monomial ideal M is played by Laurent monomials x b . 

Definition 7. (Rank of a monomial) The rank of a monomial x b with respect to M 
is one less than the minimum degree of any monomial x a that satisfies x b e (x a )\x a M. 

This definition is restricted to honest monomials x b , where b > 0. Just before the 



statement of Theorem 13 we shall extend the definition of rank to all Laurent monomials. 

The rank measures how deeply a monomial x b sits inside the ideal M. We have 
rank(x b ) > if x b e M and rank(x b ) = — 1 otherwise. Rank zero monomials form 
the border of M. Let MonSoc(M) = {x c £ M | ^x c e M \/i) denote the set of 
socle monomials of K[x]/M. See Figure [2] for a picture of a monomial ideal. The ideal 
generators are the large black circles, and monomials in M are labeled by their rank. The 
socle elements are the black squares, and other standard monomials are white squares. 

Definition 8. (Reflection invariance) A monomial ideal M is reflection-invariant if 
there exists a canonical monomial x K such that the map <p : x c h-> x K /x c defines an 
involution of the set MonSoc(M). This requires that every socle monomial divides x K . 
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Using notation as in P, §5], we note that our artinian monomial ideal M is reflection- 
invariant with canonical monomial x K if and only if the following identity holds: 

M [K+e] = (MonSoc(M)). (12) 

Here e = (1, 1, . . . , 1) and M' K+e l is the Alexander dual of M with respect to K + e. 

Definition 9. (Genus) The ideal M is level if all socle monomials have the same degree. 
If this holds then one plus that degree is called the genus of M, denoted g = genus(M). 

Example 10. Let M be the ideal generated by the seven underlined monomials in (J2]). 
Then M is level of genus g — 4, because all six socle monomials in Q are cubics, and 
M is reflection-invariant. The canonical monomial x K = x\x\x\ has degree 2g — 2 = 6. 
But, the rank of x K is equal to g — 2 = 2, as can be seen from the following lemma. □ 



For u = {ui, . . . , u rn ) G Z m we abbreviate degree" 1 " (u) = J2i- U , 



>o 



Ui. 



Lemma 11. Let M be an artinian monomial ideal. Then every monomial x b satisfies 
rank(x b ) = min degree" 1 " (b — c) — 1. (13) 

x c SMonSoc(M) 

Proof. The condition x b G (x a )\x a M in Definition [7] is equivalent to x b ~ a G K[x]\M. 
Maximizing the degree of x c = x b ~ a subject to this condition is equivalent to minimizing 
the degree of x a . But, since M is artinian, the maximal degree among its finitely many 
standard monomials is attained by one of the socle monomials x c G MonSoc(M). □ 



Remark 12. Formula (13) resembles the formula in j3j Lemma 2.2] for the rank of a 



divisor on a finite graph. We shall exploit this resemblance at the end of this section. □ 



The formula in (13) can be rewritten to be reminiscent of S-pairs for Grobner bases: 



, / lcm(x , x ) 
rank(x ) = min degree — 1. (14) 

x c eMonSoc(M) \ X c ' 



We now define the rank of an arbitrary Laurent monomial x by the formula (13). 
This is consistent with Definition [7j and it is the natural extension to monomials some 
of whose exponents are negative. The following main result gave this section its title: 

Theorem 13. Let M be a monomial ideal that is artinian, level, and reflection-invariant. 
Then M satisfies the Riemann-Roch formula, i.e., every Laurent monomial x b satisfies 

rank(x b ) - rank(x K /x b ) = degree(x b ) - genus(M) + 1. (15) 

Proof. We denote x K /x c by x c . Using the formula for rank shown in Lemma 11, the 



left hand side of (15) equals 



min degree" 1 " (b — c) — min degree" 1 " (c — b). (16) 

x c eMonSoc(M) x c eMonSoc(M) 
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Figure 2: The monomial ideal M = (x 9 ,x 6 y A ,x 5 y 7 ,x 2 y s ,y u ) is Riemann-Roch of genus 
12, with canonical monomial x 9 y 13 and Soc(M) = {x 8 y 3 , x 5 y 6 , x A y 7 , xy 10 }. Generators 
and socle monomials highlighted in dark. Monomials in M are labeled with their rank. 
The square boxes correspond to the standard monomials of K[x]/M. The dotted lines 
mark the boundary of the staircase region of M[ K+e l. Note that the identity Jl2|) holds. 
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For any socle monomial x c we have degree" 1 " (b — c) — degree" 1 " (c — b) = degree(x b ) — 
degree(x c ) = degree(x b ) — (genus(M) — 1), and hence 

degree + (b — c) = degree + (c — b) + degree(x b ) — (genus(M) — 1). (17) 



Taking the minimum of degree (b — c) over x c G MonSoc(M), equation (17) implies 
min degree" 1 " (b — c) (18) 

x c eMonSoc(M) 

= min (degree + (c — b) + degree(x b ) — (genus(M) — 1)) 

x c eMonSoc(Af) ' 

= I min degree" 1 " (c — b) J + degree(x b ) — (genus(M) — 1). 

\x c eMonSoc(M) J 



Since the map (j) is an involution, we can replace c by c in the second row of (16). It 
then follows from (18) that (16) is equal to degree(x b ) — genus(M) + 1, as desired. □ 



Remark 14. Theorem [13] can be extended to monomial ideals M that are artinian and 
reflection invariant but not necessarily level. Such M arise as initial ideals from directed 
regular (indegree = outdegree) graphs. Following [1, §2.3], we define genus min (M) as one 
minus the minimum degree of a socle monomial of M, and genus max (M) as one minus 
the maximum degree of a socle monomial of M. Using a technique similar to that in 
the proof of Theorem 13 we can derive the following Riemann-Roch inequalities: 

genus min (M) - 1 < degree(x b ) - rank(x b ) + rank(x K /x b ) < genus max (M) - 1. (19) 



Of course, the above inequality generalizes the Riemann-Roch formula ([15]): if the ideal 
M is also level then genus max (M) = genus min (M) = genus(M), and the Riemann-Roch 
formula (15) immediately follows from the inequalities ( |l~9~[ ). □ 

We say that a monomial ideal M is Riemann-Roch if it is artinian, level, and reflection 
invariant. See Figure [2] for an example in two variables. In what follows we assume that 
M is a Riemann-Roch monomial ideal. The next corollaries are formal consequences of 
the Riemann-Roch formula, as is the case for algebraic curves and graphs. 

Corollary 15. //x b is a multiple o/x K then rank(x b ) = degree(x b ) — genus(M). 

Proof. If x K divides x b , then rank(x K /x b ) = —1. Plugging this equation into the 
Riemann-Roch formula gives the assertion. □ 

Note that, by definition, the degree of the canonical monomial x K equals twice the 
socle degree. We record the following general facts about the canonical monomial x K . 

Corollary 16. The canonical monomial of a Riemann-Roch monomial ideal M satisfies 
degree(x K ) = 2 • genus(M) — 2 and rank(x K ) = genus(M) — 2. 

Experts will note that the rank is off by one when compared to the canonical divisor 



of an algebraic curve or metric graph. This discrepancy will be addressed in Remark 23 
below. We now prepare for an analogue of Clifford's theorem on special divisors. 
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Lemma 17. The rank is superadditive for monomials x a and x b of non-negative rank: 

rank(x a • x b ) > rank(x a ) + rank(x b ). (20) 

Proof. Consider an arbitrary monomial x c of degree at most rank(x a ) + rank(x b ) such 
that x c divides x a • x b . The following formulas define monomials x c and x c such that 
x c divides x a and x c divides x b : 

{a, if Ci > di, „ \ci-ai if d > di, 

and c i = < 
Ci otherwise, 10 otherwise. 

By construction, x c ' ■ x c " = x c . This implies that either degree(x c ) < rank(x a ) or 
degree(x c ) < rank(x b ). In other words, either x a /x c or x b /x c is in M, and hence, their 



product x x x is also in M. From this we infer the inequality (20) as follows: Since x c is 
an arbitrary monomial of degree less than or equal to rank(x a )+rank(x b ), and x c divides 
x a ■ x b , we know that any monomial that "defines" the rank of x a • x b (i.e., a monomial 
x d of minimum degree such that x d divides x a • x b and x& x ^ M) has degree strictly 
greater than rank(x a ) + rank(x b ). Hence, rank(x a • x b ) > rank(x a ) + rank(x b ). □ 

Corollary 18. (Clifford's Theorem) Let x b be a monomial dividing x K such that both 
rank(x b ) and rank (x K /x b ) are non-negative. Then rank(x b ) < (degree(x b ) — l)/2. 



Proof. Lemma 17 and Corollary 16 imply 

rank(x b ) + rank(x K /x b ) < rank(x K ) = genus(M) - 2. 
From the Riemann-Roch formula we have 

rank(x b ) — rank(x K /x b ) = degree(x b ) — (genus(M) — 1). 
The desired conclusion follows by adding these two identities and dividing by 2. □ 

The construction of all Riemann-Roch monomial ideals of genus g works as follows. 
We first fix a monomial x K with degree 2g — 2. Next we choose a set Ai of monomials 
of degree g — 1 that divide x K . Then there exists a unique artinian monomial ideal M 
whose socle is spanned by the monomials in M. and their complements relative to x K : 

MonSoc = M. U { x K /x b : x b G M. }. (21) 

Namely, the ideal M is the intersection of the irreducible ideals . . . , x^ +1 ) where 

x c runs over the set MonSoc. Then M is artinian, level, and reflection- invariant. 

We shall now make the connection to the Riemann-Roch theorem for graphs. As in 
Section 2, we let G denote a saturated graph on n nodes, with Uij > edges between 
nodes i and j, and M G the initial monomial ideal in K[x\, . . . , x n _i] of the toppling ideal 
Iq with respect to a reverse lexicographic term order having x n as smallest variable. 

Theorem 19. Let G be a saturated graph with n vertices, e edges, and node i having 
degree d^. Then the monomial ideal Mq is Riemann-Roch with canonical monomial 



x 



n-l 



JQ x A+u in -2 and genus ( MG ) = e _ n + 2 . (22) 
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Proof. The monomial ideal Mq is artinian, and it is level because all the socle monomials 
sj- in Corollary [5] have the same degree e — n + 1. This quantity is the cyclotomic number 
(or genus) of the graph G, which, by [5], coincides with the common degree of all maximal 
parking functions. There is a natural involution <fi on the set of (n — 1)! maximal flags 
T of subsets in [n — 1]. It takes a flag T : T\ C T 2 C ■ ■ ■ C T„_ 2 C T n _i to the reverse 
flag <p(T) : T n ^i\T n _ 2 C T rt _ 1 \T„_ 3 C • ■ ■ C T^^Ti C T n _i. Using the identification 
between flags and socle monomials in Corollary [5j we have 

S 4>(T) = x K /s r , (23) 



where x K is the monomial defined in (22). Hence, Mq is also reflection-invariant. □ 



Remark 20. Not every Riemann-Roch monomial ideal arises as an initial monomial 
ideal Mq for a connected graph G. To see this, note that the number of socle monomials 
of Mq is at most m! where m is the number of variables of Mq. On the other hand, 
Riemann-Roch monomial ideals can in general have more than m\ socle monomials. 
Furthermore, the initial monomial ideal Mq for a connected graph G is not necessarily 
Riemann-Roch. To see this, consider the four-cycle C4. Its initial ideal is Mq a = 
{x\,X2, X3) 2 , with socle monomials x\, x 2 and x%, and Mq 4 is not reflection-invariant. □ 

In the remainder of this section we show how the familiar Riemann-Roch theorem for 



graphs is derived from Theorem 19 While the proof still assumes that G is saturated, 
that hypothesis will be removed in the next section. The following algebraic definitions 
are valid for any undirected connected graph G on [n]. Here G need not be saturated. 

The Laurent polynomial ring K[x ±:L ] = Kfa^ 1 , . . . , x^ 1 } is a module over the polyno- 
mial ring K[x] = K.[x\, . . . , x n }. The Laplacian lattice image z (Ac) is a sublattice of rank 
n — 1 in Z n , and we write Lq for the corresponding lattice module, as in [9l Definition 
9.11]. Thus, Lq is the K[x]-submodule of Kfx 1 *" 1 ] generated by all Laurent monomials 
x w that have degree zero in the grading by the group Div(G). If G is saturated then 
Lq is generic and the Scarf complex is a minimal free resolution by [HJ Theorem 9.24]. 
That Scarf complex is precisely the Delaunay triangulation in [1] , and our point here is 
to redevelop the Amini-Manjunath approach in the language of commutative algebra. 

Consider the set of all Laurent monomials x c that are in the socle of the module Lq: 

MonSoc(L G ) = {x c £ Lq | X;X C G Lq Vz } . 

This socle is a set of Laurent monomials on which the lattice Lq acts with finitely many 
orbits, so the computation of MonSoc(Lc) is a finite algorithmic problem, as in 0, §9.3]. 



The problem's solution is given by the socle monomials (11 ) of our monomial ideal Mq. 



Lemma 21. The socle monomials of the lattice module Lq are precisely of the form 
sj- ■ x w /x n , where sj- G MonSoc(Mc) and x w runs over the minimal generators of Lq. 

Proof. Since G is assumed to be saturated throughout this section, the lattice module 
Lq is generic in the sense of 0, Definition 9.23], with Mq being the reverse lexicographic 
initial ideal of the corresponding lattice ideal Iq. We claim that the stated character- 
ization of the socle is valid for any generic lattice module that is artinian. Indeed, by 
the proof of [TOl Theorem 5.2], the Z"-degrees of the n-th syzygies of Lq are the vectors 
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u + w, where u runs over the Z n_1 -degrees of the (n — l)-st syzygies of Mq and w is 
any vector in the lattice. The socle degrees of Mq are the vectors u — e\ — • • — e n _i in 
Z rt_1 , and the socle degrees of Lq are the vectors u + w — e\ ■ ■ ■ — e n -\ — e n in Z n . □ 

We now identify Laurent monomials x u with divisors on the graph G. The z-th 
coordinate U{ of the exponent vector u is the multiplicity of node i in the divisor x u . 
The degree of the divisor x u is its total degree as a monomial, degree(x u ) = ui + - ■ - + u n . 



The rank of the divisor x u is defined by the same formula (13) as in Lemma 11 



rank(x u ) = min degree + (u — c) — 1. (24) 

x c eMonSoc(L G ) 

Thus, rank(x u ) > if and only if x u lies in Lq. Our definition of rank of the divisor x u 
coincides with the rank of u as in [3]. To see this, use Lemma 2.7 in [5] and note that 
the exponents of the socle Laurent monomials of Lq are the elements of the set in [3]. 

We finally define the canonical divisor of G to be the monomial 

k . . di-2 <h-2 . . . d n -2 
-X. — Xi x 2 x n , 

where di = Yuj^i u n * s the degree of node i. Finally, we recall that the genus of G is 
e — n + 1, where e is the number of edges. The following is precisely [3, Theorem 1.12]: 

Theorem 22 (Baker- Norine). Riemann-Roch holds for any divisor x u on the graph G: 

rank(x u ) - rank(x k /x u ) = degree(x u ) - genus(G) + 1. (25) 



Proof. It follows from Lemma 21 that all socle monomials of Lq have degree equal to 
the genus of G minus one. The lattice module Lq is also reflection-invariant, in the 
sense that x u G MonSoc(Lc) implies x k /x u e MonSoc(Lc). Using the representation 



in Lemma 21, the resulting involution on MonSoc(Lc) can be written as follows: 

,,d 
n 



<P((s r /x n ) ■ x w ) = (s HT) /x n ) ■ x~ w • I — ^ ) , (26) 



where 4>(T) denotes the reverse flag as in (23 ). Note that the image of is in MonSoc(Lc) 
since x d ™ j \Xl=\ X T" ^ s * n ^G- The proof of Theorem 22 is now entirely analogous to 



that of Theorem 13 In other words, our argument for the validity of the Riemann-Roch 
formula for reflection-invariant artinian level monomial ideals generalizes in a straight- 
forward manner to reflection-invariant artinian level lattice modules. □ 

Remark 23. The rank of the canonical monomial of Mq equals the rank of the canonical 
divisor of the graph G, but the degree of the former is two more than that of the 
latter. □ 



4 Non-Saturated Graphs 

We turn to graphs G that are not necessarily saturated. The binomials in ([6]), their 
syzygies in (j9), and the ideals Iq and Ma are still well-defined. However, the minimality 
in Theorem [2 is no longer true, and the choice of term order is more subtle, as the next 
example shows. 
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Example 24. Let G be the edge graph of a triangular prism, labeled so that Iq equals 



[(a 3 - bed, b 3 - ace, c 3 - abf, d 3 - aef, e 3 - bdf, f - cde) : (abode f)°°]. 

This toppling ideal has 22 minimal generators, and its free resolution has Betti numbers 
(1, 22, 92, 147, 102, 26). The same holds for the ideal that represents parking functions: 

M G = in(J G ) = (a 3 ,abc,ace 2 ,b 3 ,b 2 e 2 ,c 3 ,cde,d 3 ,d 2 e 2 ,e 3 a 2 b 2 ,a 2 c 2 ,a 2 d 2 , 

abde, b 2 c 2 , bed 2 , a 2 be 2 , a 2 de 2 , ab 2 d 2 , ac 2 d 2 , b 2 d 2 e, bc 2 e 2 ). 

Here the reverse lexicographic term order with a>b>c>d>e>f was used. However, if 
we take e>d>c>b>a> f then the reverse lexicographic Grobner basis requires two more 
binomials. Now, the initial monomial ideal has Betti numbers (1, 24, 98, 153, 104, 26). □ 

To explain the phenomenon in this example, we fix a spanning tree T of the graph G 
that is rooted at the node n, and we order the unknowns according to a linear extension 
of T. Thus, we fix an ordering of [n] such that i > j if the node i is a descendant of the 
node j in T. A term order on K[x] is a spanning tree order if it is a reverse lexicographic 
term order whose variable ordering is compatible with some spanning tree rooted at n. 
One spanning tree order is the toppling order considered in Theorem 10]. See also 
[TT| §5] for a discussion of Grobner bases of toppling ideals in the inhomogeneous case. 

Theorem 25. The toppling ideal Iq is generated by the binomials x /_i>J — x J ^ 7 where 
(I, J) runs over splits of [n] such that the subgraphs of G induced on I and J are con- 
nected. For any spanning tree order, these binomials form a Grobner basis of Iq with 
initial monomial ideal Mq. The complexes constructed in Isl) are cellular free resolutions. 



Proof. The first paragraph in the proof of Theorem [2] is valid in the non-saturated case. 
It shows that the binomials x 1 ^ 7 — x J_i>/ lie in Iq. For the spanning tree term order, the 
leading monomials are x /_i>J , where n G J\I, and hence the initial ideal in(/c) contains 



the monomial ideal Mq of (10). Again, both ideals are artinian of the same colength 
in K[x\ n ], and hence they are equal. This establishes the Grobner basis property. The 
argument in the proof of [HI Theorem 14] shows that the property that the subgraphs of 
G induced on I and J are connected characterizes a minimal Grobner basis of Iq and 
the minimal generators of Mq. In particular, these binomials x /_ ^ J — x" 7 ^ generate Iq. 

Our last assertion states that (|8f with differentials ^ gives a free resolution of Iq, 
and dropping the last term in ^ gives a free resolution of Mq. This claim is proved 
by deformation to generic monomial modules, as explained in P, §6.2]. To be precise, 
in our situation we replace G by a nearby saturated graph G e with fractional edge 
numbers Uij(e) between any pair of nodes. The monomial ideal Mq e is generic and 
degenerates to M G . The lattice ideals Iq and Iq c are represented by the corresponding 
lattice modules Lq and Lq c . These are submodules of the Laurent polynomial ring as in 
P Definition 9.11]. The lattice module Lq c is generic and degenerates to Lq. According 
to [HI Theorem 6.24], the Scarf complex of Mq € with labels from G gives a free resolution 
of Mq. Likewise, the Scarf complex of the generic lattice module Lq c with labels from G 
gives a free resolution of Lq. The resulting minimal free resolution of Iq c degenerates to a 
(typically non-minimal) resolution of Iq, using [21 Corollary 9.18]. These free resolutions 
of Mq and Iq are cellular because they are given by labeled simplicial complexes. □ 
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Here is an example that illustrates the degeneration used in the proof above. 



Example 26. Let n = 4 and G be the 4-cycle 1 — 2 — 3 — 4 — 1. For 5, e G N consider 
the graph Gs, e that has 5 edges for every edge in G and e edges for every non-edge of G. 
Then Gg, e is saturated for 5, e > 0. The Scarf complex in Example [3] gives the minimal 
free resolution of Mq Sc = in(J(j 5e ) and this lifts to the minimal free resolution of Ig 6c - 
By Theorem 19, the monomial ideal Ma St is Riemann-Roch, its genus is 4<5 + 2e — 2, and 
its canonical monomial equals x K = xl 6+e ~ 2 x 2S+2t ~ 2 xl s+e ~ 2 . The involution x b t— y x K b 
on its six socle monomials is given by swapping the two rows below: 

( x 5 ~ 1 x 5+€ ~ 1 x 25+e ~ 1 x 5 ~ 1 x 2S+e ~ 1 x 5+e ~ 1 x 5+e ~ 1 x 25+e ~ 1 x 5 ~ 1> \ 

MonSoc(M G ,J = k-lJ+el>l' JJ+e-Le-L %S-1 ' Vl el Ji+el >• (27) 
l^l x 2 x 3 ) x l x 2 x 3 5 x l x 2 x 3 J 

Setting 5 = 1 and e = 0, we get the parking function monomial ideal of the 4-cycle 

M Gl0 = Mq = (x 1 ,x 2 ,x 3 ) 2 = (x 2 1 ,x 1 x 2 ,x 1 x 3 ,xl,x 2 x 3 ,xl). 

Here, MonSoc(Mc) = {xi,x 2 ,x 3 }. This ideal is not reflection-invariant and hence not 
Riemann-Roch. The cellular resolution of Mq induced from Mq s e is not minimal. □ 

We now take a closer look at the combinatorial structure of our resolutions. Let 
Bary(G) denote the first barycentric subdivision of the (n — 2)-simplex, whose 2 n ~ 1 — 1 
vertices, namely the non-empty subsets I of [n — 1], are labeled by the corresponding 
monomials y^^wv. Thus, Bary(G) is the cellular free resolution of Mq = m(Ia) referred 



to in Theorem 25 Each simplex in Bary(G) is labeled by the least common multiple 
of the monomials that label its vertices. For any c 6 N n_1 we write Bary(G)^ c for the 
subcomplex consisting of all simplices in Bary(G) whose labels properly divide x c . 

Corollary 27. The number of minimal i-th syzygies of the monomial ideal Mq in degree 
c is equal to the rank of the reduced homology group ifj_i(Bary(G)^ c ; IK). 



Proof. This follows immediately from Theorem 25 and 0, Theorem 4.7]. □ 



We next state the analogous result for the lattice module La, that is, the K[x]-module 
generated by all Laurent monomials whose exponent vector lies in the Laplacian lattice 
image z (Ac). We identify this lattice with Z n /Ze by writing its elements as Aq ■ v where 
each v e Z™ is unique modulo Ze = ker(Ac). The tropical metric on Z n /Ze is 

dist(u, v) = maxj \ui + Vj — Uj — Vi\ : 1 <i < j < rt}. 

We write Apt(G) for the corresponding flag simplicial complex. Thus, Apt(G) is the 
simplicial complex whose simplices are subsets S of Z n /Ze such that dist(u, v) < 1 for 
u, v G S. The notation "Apt" refers to the fact that this infinite simplicial complex is 
the standard apartment in the affine building of Lie type A n -\. It is well-known that 
Apt(Cr) is pure of dimension n — 1 and that it triangulates the (n — l)-dimensional affine 
space ]R n /]Re. For more on buildings and their connection to tropical geometry, see [7]. 

The apartment Apt(G) is precisely the same as the Delaunay triangulation con- 
structed in P], and it also coincides with the Scarf complex of G e that we used to 
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give a cellular resolution of Lq. The number of i-faces of Apt(G) modulo the lattice 
action is given by ([7]). Each vertex v of Apt(G) is labeled by the corresponding Laurent 
monomial x Agv , and each face is labeled by the least common multiple of its vertex 
labels. Thus, each face of Apt(G) is labeled by a Laurent monomial of degree > 0. We 
write Apt(G)_< c for the subcomplex of all simplices whose label properly divides x c . 

Corollary 28. The number of minimal (i + l)-st syzygies of the lattice module Lq in 
degree c is the rank of the reduced homology ifj(Apt(G)_< c ; IK). The sum of these ranks 
over all c modulo image z (Ac) counts the minimal i-th syzygies of the toppling ideal Iq. 

We conjecture that the ranks of the homology groups in the two corollaries coincide. 

Conjecture 29. The Betti numbers of the toppling ideal Iq do not increase when passing 
to the initial ideal M G . More precisely, for all i > and all c G N n_1 ; we have 

^_ 1 (Bary(G)^ c ;K) ~ ^(Apt(G)^ c ; K). (28) 
This conjecture has been verified for many graphs using the software Macaulay2. We 



A 



note that the two simplicial complexes appearing in (28) are different from the complex 
Ap used in Hochster's formula for the Betti numbers of a lattice ideal [TT| Theorem 7.4]. 

Example 30. The simplicial complexes Apt(G)^ c can be large even for small graphs. 
Let G be the graph on four nodes, labeled a, b, c, d, with Laplacian matrix 

/M12+U13+M14 -U 12 ~«13 ~«14 \ 

—Ui2 U12+U23+U24 -«23 -U 2 4 
-Ui3 -U 23 U 13 +U 2 3+U U -U34 

\ -U U -U 2 4 -M34 Ul4+U 2 4+U 34 J 

Both the toppling ideal and the ideal of parking functions are complete intersections: 

I G = ( a 2 -b 2 ,b- c,c 3 - d 3 } and M G = (a 2 ,6,c 3 ). 

The monomial ideal Mq has one minimal first syzygy in degree c = (2,0,3,0). The 
simplicial complex Bary(G)^ c consists of two isolated nodes a 2 and c 3 . The simplicial 
complex Apt(G)^ c is two-dimensional but it has the homology of a circle. It consists 
of 12 triangles, 28 edges and 16 vertices, labeled by the following generators of Lq\ 



( 


2 


-2 





o\ 




-2 


3 


-1 










-1 


4 


-3 


\ 








-3 


3/ 



c c 2 c 3 



a 2 c a 2 



a 2 c 2 



a 2 a 2 c 3 



a 2 c a 2 c 2 a 2 c 3 



' 6' 6 2 ' 6 3 ' a 2 ' a 2 6' d 3 ' d 3 ' c 2 ' bd 3 ' be b 2 d 3 ' b 2 ' b 3 ' b 4 ' b 5 ' 
The lattice module Lq has one second syzygy in this degree, translating into a first 
syzygy of /g- It is represented in Apt(G)^ c by the 4-cycle l,a 2 /b 2 ,a 2 c 3 /b 2 d 3 ,c 3 /d 3 . □ 

At present, no explicit minimal free resolution of Mq is known. Finding such a 
resolution was stated as an open problem by Postnikov and Shapiro in [12, §6]. We do 
not even know whether the Betti numbers depend on the characteristic of the field K. 

An explicit formula for the Betti numbers of the toppling ideal Iq was conjectured 
by Wilmes in [15]. See also (TTJ §7.4]. Wilmes' formula is combinatorial, and it has 
been verified for all graphs with n < 6 nodes. At present we do not know how to relate 



Wilmes' conjecture to the ranks of the homology groups in Corollaries 27 and 28 



It is known, thanks to [T3J Theorem 3.10], that Conjecture 29 is true for the maximal 
syzygies, with index % = n — 2. We have the following combinatorial characterization: 
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Corollary 31. The maximal syzygies of the parking function ideal Mq, or of the toppling 
ideal Iq, are in bisection with the acyclic orientations of G with node n as unique sink. 



See [TT| Theorem 7.6] for an alternative but equivalent formulation of this result. 

Proof. The monomial ideal Mq is artinian, so its maximal syzygies correspond to the 
socle elements. These are the maximal parking functions, and, by [5j Theorem 4.1], they 
correspond to cyclic orientations of G with node n as unique sink. Since all maximal 
syzygies of M G = in(/ G ) lift to maximal syzygies of Iq, the same result holds for I G . □ 

We now derive the Riemann-Roch theorem for non-saturated graphs G. Let Mq be 
the initial ideal with respect to a spanning tree order on the variables with x n as the least. 
By Corollary |3~lj we know that the socle monomials of the Laurent monomial module Lq 
are s ■ x^ 1 ■ x w where s runs over all socle monomials of Mq and x w runs over minimal 
generators of Lq. Unlike in the saturated case, the monomial ideal Mq is generally 
not reflection- invariant. But the Laurent monomial module Lq is always reflection- 



invariant. To see this, we use Lemma 3.2 of [3] to deduce that sq-jx n , defined in (11), is 
not contained in Lq for any complete flag T of [n]. This implies that sq-/x n is a socle 
element of Lq, since every Laurent monomial of degree greater than degree (sf/x n ) = 
genus(G) — 1 is in Lq- We now immediately verify that Lq is reflection-invariant with 
the involution on MonSoc(LcO that takes srx w /a; n to s^-r^isr™ x^~ l / YYiZi x T n > where 
</>(T) is the reverse flag of T exactly as in the proof of Theorem [22] The canonical 
monomial is 

Y k _ Ji-2 Jj-2 ™dn-2 
A. x 2 " " " Jj n , 

where di = Uij is the degree of node i. Hence, Lq satisfies the Riemann-Roch 



formula, in its monomial formulation (15), with M = Lq and genus(M) = genus(G). 
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